Optimal Control of One-Qubit Gates 
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Abstract 

We consider the problem of carrying an initial Bloch vector to a final Bloch vector in a specified amount of 
time under the action of three control fields (a vector control field) . We show that this control problem is solvable 
and therefore it is possible to optimize the control. We choose the physically motivated criteria of minimum 
energy spent in the control, minimum magnitude of the rate of change of the control and a combination of both. 
We find exact analytical solutions. 



>1 Introduction 
O 

^^ecent advances in experimental physics allowing for manipulation and measurement of single quantum systems, 
"iiave stimulated a flurry of investigations on the control of quantum systems, and more or less formal schemes have 
^—ieen advanced[l, 2[ 3, ^. The conditions under which a given quantum system is completely controllable have 
^^een explored js, ^, 7| and some limits of quantum controllability j^] have been found. Quantum control theory 
fTiave several important applications including quantum state engineering ||l0|] , control of chemical reactions [pt [l^ , 
"^aser cooling of molecular degrees of freedom quantum register initialization jl^ and the fabrication of 

CTcj ibust quantum memories [|l6[ . 

^ ' , , A major application of the theory of quantum control is the subject of quantum computation. The physical 
implementation of a quantum computer is a major challenge and many proposals|^^ including ion traps, optical 
C^8,vities, and quantum dots have been advanced. Promissing practical implementations should be scalable, facing 
gjfie problem of heat dissipation which gets worse along with the shrinking of the size of the proposed physical 
. system. In this work we address the problem of carrying an initial qubit (more precisely of an initial Bloch vector) 
. specified final qubit in a given amount of time using the minimum amount of energy possible. This paper 

^^s, organized as follows: first we review the controllability of Bloch vector, and then we formulate and solve the 
V^jroblem of optimal control under the criteria of minimum energy, minimum energy derivative and a combination 
. 5^fboth. 

Two level systems adequately model many physical systems (spin 1/2, photon polarization, atoms in (quase) 
monochromatic electromagnetic fields, etc), despite its simplicity. A general 1-qubit gate can be represented as a 
2-level quantum system, and the most general Hamiltonian for such a system can be written as H = hoi -\- h • a, 
where a are Pauli's matrices, /iq determines the zero energy reference, and /i is a classical vector. Since we can 
use any h we want, and this is the field we use to control our system, here on we call h our vector control. We 
moreover assume = Q. 



2 Bloch vector Controllability 

All of the information of the quantum state of a two level system is completely determined by its density matrix 
p, or equivalently by its Bloch vector s{t) = ^ Ti{pd). The dynamics of Bloch vector, given by well known Bloch 

equation s(t) = b x s, where b = 2h/h, can be put in the more explicit way s = {b^Jx + byjy + bzJz)s, where the 
Js are the rotation generators. In this case we can formulate the problem of taking an initial state pi to a final 
state in a specified amount of time T = tf — ti. After rescaling, we take ti = 0,tf = 1. Remember that we 
assume that all three components of h (or of b) are control fields. Thus, since we have all three rotation generators, 
the system is completely controllable in the sense that every rotation can be reached from the identity |@, ||] . In 
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other words, any final vector can be reached from any other initial vector (of the same length, or the same degree 
of mixture) in any finite time, provided there are no constraints on the size of the control fields. Moreover, the 
motion equation can be inverted [p^ to give b — sx ss~^ — f{t)s, where fit) is an arbitrary function, which shows 
that not only the control problem is solvable, but that it is solvable even when a path to realize such control is 
given. Moreover, even in this event the solution is not unique. 



3 Optimal Control 

As mentioned before, the larger the fields used for control the greater the amount of heat to be dissipated. Then, 
proved the complete controllability of the one qubit gate, it is meaningful asking which is the control vector field 
required to perform an arbitrary rotation operation with the minimum amount of energy spent. We address this 
problem an optimal "classical" control problem for the Bloch vector, that is we shall extremize the cost functional 

S = f ' dtllb»b + X» (k-bx s)\ , (1) 



where b = 2h/h, and A is a (vector) Lagrange multiplier, and the momentum corresponding to s. Notice that the 
form of the cost functional above is not arbitrary: were the vector control field a magnetic field acting on a spin 
half system, or an electromagnetic field for a charged two level system, the (electro) magnetic energy would have 
had the assumed form. For the sake of simplicity we scale the variables to get the following cost functional 

S = dtl^^b»b + X» (k-bx s^^ , (2) 

where the constant a have been introduced to have s and A in the same units. 
Not all of the resulting Euler-Lagrange equations 

b — asxX, s = bxs, \ = bxX, (3) 

are true dynamical equations: the first equation is a constraint. A few simple calculations show that the system of 
equations (||) possesses the following constants of motion = s • s, A^ = A • A and i/ = X • s = Ascos(6'), where 9 
is the angle between A and s. The equations of motion can be put in the equivalent form 




cos(6') 1 \ 



whose solution reads 



d[at]\Xj \ -X^ X-s \Xj \ -1 cos{9)J\X 



_ siii{9 - at sin 9) sin(aisin0) ^ /^^('^)^ 

x) ~ sm{9) -sin(aisin6l) sin(6i + at sin 6*) ) \X{0)) 



(4) 



(5) 



where 



We have assumed that s and A are unitary vectors. Notice that we have obtained an orthogonal transformation. 
In fact this could have been anticipated by showing that b is constant both in norm and direction. Let us see that 
the solution is valid for any value oi 9, so we have freedom to choose 9 to our best convenience. For example, we 
can choose A(0) = s(l). It follows cos 61 — s(0) • s(l), and sin(a sin(6')) = sin(6'). This choice leads to the solution 

sm{b{l-t)) ^ smjbt) ^ 
sm(6) sm(6) 

7', ^ N ArcSin (sin(^)) ^ , ^ ^ , , ^ 

= KO) = -±^Jlg^ X Sf ^{9 + 2^n)si = bi{t; n)s^. (7) 

sm^f j 

Equations (||) and (|^) constitute the solution to the proposed problem, and agree with the result of reference JTof . 
Let us see that in the scaled variables the magnitude of the field is essentially equal to the value of the angle 
between the initial and the final Bloch vectors. The solution corresponds to a rotation around {—)s± at constant 
speed, from the initial Bloch vector to the final one. We see then that the solution is not unique: due to the 
multivaluedness of the function ArcSin, there is an infinite number of solutions. Each solution corresponds to a 
local minimum, the global minimum being the solution whose magnitude is equal to the angle between the initial 
and the final Bloch vector. The other solutions correspond to an integer number of turns followed by this angle. 
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Finally, the other solutions correspond to reaching the final Bloch vector from the initial one but rotating the other 
way. Choices different to A(0) — s(l) are also meaningful, but should lie on the plane which contains the initial and 
final Bloch vectors. Let us see that the solution found breaks down when the initial and final vector are antipodal. 
In this case there are an infinite number of solutions all of them spending the same amount of energy. In this case 
the choice of an initial Lagrange multiplier vector which is not (anti)parallel to the initial Bloch vector leads to a 
particular solution. 

Now, we remark that the minimum energy control corresponds to minimum path on the Bloch sphere. The 
length I transversed by the tip of the Bloch vector is given by 

1= T '^^= bxsdt^s b{t)\ ain{e {t))\dt < s / b{t)dt. (8) 
Jo Jq Jq Jq 

This means that, for a given magnitude of the control field, the transversed length is maximized when the control 
vector field makes a right angle with the Bloch vector. In other words, if we restrict to control fields that are 
perpendicular to Bloch vector, and if Bloch vector transverses a longer path, it is necessary to have a larger average 
control field. This remark allows us to restrict ourselves to consider the vector control to point out in the direction 
of Si X Sf . Making the Ansatz 

m = Kt)s±, (9) 

after some algebra we obtain 

sin(g-/>ft'K) , , Mj^b{t')dt') ^ 
= ^) + sin(^) (1") 

were b should satisfy the equality b{t)dt = 9. If we write b{t) as 9 + 5{t) we see that the average value of 5{t) over 

the unitary interval is zero and that /p 9{t)dt = 0"^ + S'^(t)dt. Thus, the minimum is attained when the control 
vector field is constant from the initial till the final time. It is possible to give an alternative argument which 
shows that the solution of minimum fluence is the same as the shortest geodesic. Inverting the Bloch equation we 

obtain the control field b — sx s + fs, and the energy spent in the control b^{t)dt = (i 's)^ ^ (s*- s)^ + /^^ dt. 
Since the second subintegral term is identically zero, and the function f{t) should be zero for the extrema, we see 
that the fluence minimization and geodesic minimization (see Eq. |^) are almost the same, and reach their extrema 
together. Had we chosen the squared length instead of the length, both expressions would have been identical. 

Since the times to perform quantum computation are generally short in low dimensional condensed matter 
systems, which are the most promising candidates, one also should analize possible limitations set by the rate at 
which control fields can be set. In particular, it is worth paying attention that the solution is a discontinuous 
one, zero before the initial time, constant between the initial and final times, and zero again from the final time on. 
Had we used the square of the time derivative of the control field instead of the square of the field itself, defining 
the cost functional 



the solutions obtained above would have been also solutions of the new problem. In this case a whole set of new 
solutions arise, which are of constant magnitude but whose direction changes with time. It is easy to construct 
such a kind of solutions. Let {b{t),s{t)} a solution of the new problem, but with a time dependent b'^{i), then 

{b{t) = b{t) + f{t)s,s}, is also a solution, no matter how the function / is chosen. In particular, we can adjust 
/ to obtain b^ a constant. For instance, if we set s{t) = cos{<j>{t))so{t) + sm{(j>{t))s±, where sq is the solution 
for the problem of minimum fluence, and cos((j)(t)) a function with value 1 both at < and at t = 1, we have 
the control field b{t) = 9 cos'^{(j))s± — (pSr — 6* sin(0) cos(0)so where Sr — s± x sq is a unitary vector needed to 
define a time dependent right triade {sq, Sr, s±}. One can choose f as f — ±\J B'^ — ff^ cos^((/)) — with the 
maximum value of so at instants where the maximum is attained, / vanishes. For the sake of definiteness we use 
(^{t) = QyttiX — t) which yields = 6^ {sin^ {(j){t)) + /i^t(2 — t)), and produces a new constant norm vector control 

b with magnitude + /i^. This solution, of course, spends more energy than the found before to perform the 
same control. 

Solutions with vanishing magnitude at the initial and final instants of time also exist. Observe that in this case 
the equation for b is 

b=-n'^sx\ (12) 
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which is a true dynamical equation. This allows for some extra flexibility: now we can add initial and final 
conditions on the value of the control field. From the point of view of the energy injected to the system, the most 
physically sensible conditions are those of vanishing control field both at the initial and the final times. We observe 
that in this case the solution should follow the shortest geodesic between the initial and the final Bloch vectors. 
In fact, since the control field begins and ends with a vanishing value, it should grow and decrease as slowly as 
possible but fast enough as to reach the maximum value necessary to have an average magnitude of at least 0. 
If it had grown to a greater value it would have grown at a larger pace, so it would have been not the minimum 
solution sought. 

Differentiating the equation (12) we obtain ^ = bx which immediately tells us that the second derivative 

of the vector control has constant norm. The first integral of this third order equation, ^ — ^ ^ ^ ~ fl^Si x s/, 
where we have assumed A(0) — Sf, can be solved under the assumption that the control vector is a second degree 
polynomial in t, with the result b{t) = £LL5HlWit(l — i)s^ = b2{t)sj_ . Numerical solution of these equations, 
without the ansatz made above also lead to the same solutions. Before expressing b completely in terms of Si and 
s/, we proceed to discuss the more general physical criterion in which one is interested on energy saving but with 
a limited rate of change of the vector control, through the cost functional 

1 db db\ r r a] 

The experience gained with the previous examples shows that the solution control tied should point (anti)parallel 
to s± . Some algebra leads to the solution 

. ^, / coshtj(t - i)\ 
m=a\MO)\[l~ cosh(l) j-^^^^--^^)--^- 

We notice that all of the solutions so far found have the form of equation (H) and therefore have the solution (|l0|). 
We only have to take care of the final value of s. This leads to the following more explicit forms for b{t) 

M«;n)^6(. + 2,„),(i-,). w„„.^l^(i-:=5i^)^ (14) 

For the second case considered, the intuitive choice, b = 7r0sin(7rt)/2 produces a value of the cost functional only 
1.5% above that of the optimal solution. Finally, for control fields of the form of equation the cost functional 
can be written in purely geometric terms. If we set = b{t), 0(0) = 0, then S can be expressed as 




where b is the average magnitude of the control field, and the accumulated angle (or the arc length) transversed 
by the Bloch vector. Equation (15), just like equation (p^), contains the other two cases: the first in the limit 



1/lu ^ 0, and the second in the limit 1/a — > but with auj'^ — fl^ fixed. Of course, bi(t; n) = lmii/^_,Q b3{t] n) and 
b2{t;n) = limi/o^o,ac:;2^n2 h{t\n). 

We have formulated and solved in an analytic way, the problem of rotation of the Bloch vector (which charac- 
terizes completely the state of a two level system) from a prescribed initial vector to a prescribed final vector, in a 
given amount of time, using an optimal control scheme which minimizes the energy spent by the control fields, or 
the magnitude of the rate of change of the control fields or a linear combination of both. We have found control 
fields perpendicular to both the initial and final Bloch vectors, and multiple local minima corresponding to arrival 
from the initial to the final Bloch vector in one or other senses or after one or more complete turns. 

This work was partly funded by DIB-UN (Division de Investigaciones, Sede Bogota, Universidad National). 
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